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INTRODUCTION 

The  goals  of  the  first  year  of  this  training  grant  were  to  develop  the  foundations  for  a  new  medical 
imaging  modality,  now  called  Radioluminescence  Tomography  (RLT),  which  utilizes  biologically 
compatible  phosphor  nanoparticles  to  produce  images  of  molecular  breast  cancer  contrasts.  The 
development  of  this  modality  is  significant  to  the  early  detection  of  breast  cancer,  because  it  provides 
additional  disease- specific  molecular  information  to  X-ray  mammography,  which  is  the  most 
widespread  medical  imaging  modality  for  breast  cancer  screening.  In  addition,  this  modality  has  the 
potential  to  provide  feedback  to  guide  Radiation  Therapy,  for  such  breast  cancer  applications  as 
Intraoperative  Radiation  Therapy.  The  first-year  research  goals  were  to  develop  a  mathematical 
framework  and  the  instrumentation  to  enable  this  modality,  and  to  fabricate  and  fully  characterize 
phosphor  nanoparticles  for  RLT.  In  addition,  this  research  grant  has  aided  in  the  development  of  a  new 
more  biologically  compatible  radioluminescent  phosphor.  Barium  Yttrium  Fluoride,  which  has  been 
thoroughly  characterized,  and  successfully  imaged  in  phantoms  and  a  mouse.  In  addition,  this  grant  has 
provided  the  funding  to  devise  a  more  sophisticated  version  of  RLT,  Limited- Angle  X-ray 
Luminescence  Tomography,  which  has  the  potential  to  achieve  high  resolution  with  a  decrease  in  X-ray 
dose  to  healthy  tissue.  For  the  first-year  training  goals,  this  grant  has  provided  for  extensive  study  in 
molecular  targeting  of  cancer,  molecular  imaging  modalities,  and  opportunities  to  engage  physicians. 

As  of  today,  the  results  of  this  grant  are:  5  journal  publications,  including  2  first  author,  with  2  more  in 
preparation;  8  conference  abstracts,  including  3  first  author,  and  2  conference  presentations;  4  courses 
taken,  including  BioE222:  Molecular  Imaging,  Med374:  Medical  Device  Design  (taken  after  the  grant 
was  accepted,  but  before  it  was  funded),  the  Comprehensive  Cancer  Research  Training  Program  at 
Stanford,  and  SIE:  The  Stanford  Institute  for  Entrepreneurs,  and  mentoring  of  1  undergraduate  and  1 
high  school  student.  The  future  for  this  grant  looks  bright,  as  this  technique  will  soon  be  tested  in  small- 
animal  models  in  vivo ;  the  next,  more  clinically  feasible  interventional  version  is  under  development. 

BODY 

Research  Accomplishments:  SOW  Aim  1:  Nanophosphor  characterization 
The  feasibility  of  X-ray  luminescence  with  Gadolinium  oxysulfide  particles  was  investigated.  We 
examined  the  practical  aspects  of  this  new  modality,  including  phosphor  concentration,  light  emission 
linearity,  detector  damage,  and  spectral  emission  characteristics.  Finally,  the  contrast  produced  by  these 
phosphors  was  compared  to  that  of  X-ray  fluoroscopy. 

Gadolinium  and  lanthanum  oxysulfide  phosphors  doped  with  terbium  (green)  or  europium  (red) 
were  studied.  The  light  emission  was  imaged  in  a  clinical  X-ray  scanner  with  a  cooled  CCD  camera  and 
a  spectrophotometer;  dose  measurements  were  determined  with  a  calibrated  dosimeter.  Using  these 
properties  in  addition  to  luminescence  efficiency  values  found  in  the  literature  for  a  similar  phosphor, 
minimum  concentration  calculations  were  performed.  Finally,  a  2.5cm  agar  phantom  with  a  1cm- 
diameter  cylindrical  phosphor- filled  inclusion  (diluted  at  lOmg/ml)  was  imaged  to  compare  X-ray 
luminescence  contrast  with  X-ray  fluoroscopic  contrast  at  a  superficial  location. 

Dose  to  the  CCD  camera  in  the  chosen  imaging  geometry  was  measured  at  less  than 
0.02cGy/sec.  Emitted  light  was  found  to  be  linear  with  dose  (R  =  1)  and  concentration  (R  =1). 
Phosphor  emission  peaks  were  less  than  3nm  full- width  at  half-maximum,  as  was  expected  from 
lanthanide  dopants.  The  minimum  practical  concentration  necessary  to  detect  luminescent  phosphors 
was  dependent  on  dose;  it  was  estimated  that  sub-picomolar  concentrations  are  detectable  at  the  surface 
of  the  tissue  with  typical  mammographic  doses,  with  the  minimum  detectable  concentration  increasing 
with  depth  and  decreasing  with  dose.  In  a  reflection  geometry,  X-ray  luminescence  had  nearly  a  430- 
fold  greater  contrast  to  background  than  X-ray  fluoroscopy. 
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The  outcome  of  this  aim  was  a  successful  demonstration  and  feasibility  assessment  of  this 
modality  using  Gadolinium  Oxysulfide.  However,  due  to  difficulties  in  making  these  phosphors 
biocompatible,  the  phosphors  were  switched  to  BaYF4,  which  may  be  more  easily  bioconjugated  and 
tagged  to  molecular  markers.  Indeed,  in  a  preliminary  study,  these  phosphor  were  targeted  to  the  Folate 
receptor  (commonly  expressed  in  breast  cancer),  and  uptaken  by  live  cells  (data  not  shown). 

Training  Accomplishments:  SOW  Aim  1:  Nanophosphor  characterization 

In  performing  these  tasks,  the  investigator  has  been  exposed  to  the  field  of  molecular  imaging,  a  new 
direction  for  this  PI.  This  research  education  has  been  aided  with  participation  in  BioE222:  Molecular 
Imaging,  which  brought  together  the  top  molecular  imaging  faculty  at  Stanford  to  teach  aspects  in  the 
hardware,  chemistry,  and  biology  of  molecular  imaging.  In  addition,  the  PI  was  exposed  to  nanoparticle 
fabrication,  including  the  processes  in  making  nanoparticles  stable  in  human  serum  with  low  toxicity. 
Also,  the  PI  gained  knowledge  in  molecular  targets,  and  the  advantages  and  disadvantages  of  targeting 
to  peptides,  hormones,  antibidies,  affabodies,  and  other  targeting  agents.  This  program  has  been  aided 
by  working  at  benchside  with  these  materials  scientists,  biologists,  and  nuclear  imaging  experts. 


Research  Accomplishments:  SOW  Aim  2:  Imaging  Hardware  Development 

This  aim  extended  the  ability  of  Radioluminescent  Imaging  with  the  construction  of  a  fully  functioning 
system,  shown  in  Figure  1,  which  has  full  spectral  and  imaging  capabilities  to  image  small- animals 
injected  with  nanophosphors  with  high  resolution.  In  addition,  this  aim  resulted  in  the  development  of  a 
novel  algorithm  which  incorporated  a  diffuse  optical  photon  propagation  model  into  the  reconstruction 
algorithm  to  recover  unresolved  dimensions  in  an  X-ray  limited  angle  (LA)  geometry.  This 
sophisticated  technique  enables  such  applications  as  image-guided  surgery,  where  the  ability  to  resolve 
lesions  at  depths  of  several  centimetres,  which  can  be  the  key  to  successful  resection.  In  addition,  the 


Figure  1:  Small- animal  imaging  box 
fabricated  to  provide  proof-of-concept  and  aid 
in  technique  and  nanoparticle  development. 


internal  organs  need  not  be  irradiated.  The  hybrid 
X-ray  /  diffuse  optical  model  was  formulated  and 
demonstrated  in  a  breast-sized  phantom,  simulating 
a  breast  lumpectomy  geometry.  Both  numerical 
and  experimental  phantoms  were  tested,  with 
lesion-simulating  objects  of  various  sizes  and 
depths.  Results  showed  localization  accuracy  with 
median  error  of  2.2mm,  or  4%  of  object  depth,  for 
small  2-14mm-diameter  lesions  positioned  from 
lcm  to  4.5cm  in  depth.  This  compares  favorably 
with  fluorescence  optical  imaging,  which  is  not 
able  to  resolve  such  small  objects  at  this  depth. 

The  recovered  lesion  size  had  lower  size-bias  in  the 
X-ray  excitation  direction  than  the  optical 
direction,  which  was  expected  due  to  the  optical 
scatter.  However,  the  technique  was  shown  to  be 
quite  invariant  in  recovered  size  with  respect  to 
depth,  as  the  standard  deviation  was  less  than 
2.5mm.  Sensitivity  was  a  function  of  dose; 
radiological  doses  were  found  to  provide  sufficient 
recovery  for  pg/ml  concentrations,  while  therapy 
dosages  provided  recovery  for  ng/ml 
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concentrations.  Experimental  phantom  results  agreed  closely  with  the  numerical  results,  with  positional 
errors  recovered  within  8.6%  of  the  effective  depth  for  a  5mm  object,  and  within  5.2%  of  the  depth  for  a 
10mm  object.  Object  size  median  error  was  within  2.3%  and  2%  for  the  5mm  and  10mm  objects, 
respectively.  For  shallow-to-medium  depth  applications  where  optical  and  radio-emission  imaging 
modalities  are  not  ideal,  such  as  in  intra-operative  procedures,  this  new  technique,  LAXLT,  may  be  a 
useful  tool  to  detect  molecular  signatures  of  disease.  For  this  aim,  the  software  and  hardware  proof-of- 
principle  experiments  were  performed. 

Training  Accomplishments:  SOW  Aim  2:  Imaging  Hardware  Development 
In  performing  these  tasks,  the  investigator  explored  various  methods  of  multimodality  imaging,  by 
fusing  X-ray  and  optical  photon  propagation  models.  This  work  was  carried  out  under  guidance  from 
Dr.  Fei  Xing,  as  proposed,  who  is  an  expert  in  X-ray  modeling  in  tissue.  A  numerical  photon 
propagation  model  was  built,  and  integrated  into  the  Pi’s  existing  code  for  optical  photon  modeling. 

The  system  development  for  pre-clinical  imaging  was  a  first  for  this  PI:  a  small-animal  system  devoted 
to  testing  modality  feasibility,  and  one  that  may  have  potential  for  high-resolution  small-animal 
imaging. 


KEY  RESEARCH  ACCOMPLISHMENTS 


•  Established  the  feasibility  of  Radioluminescent  Tomography  in  reflection  mode. 

•  Developed  a  hybrid  X-ray /Optical  photon  propagation  model  to  perform  high  resolution 
Radioluminescent  Tomography,  in  an  appropriate  geometry  for  Intraoperative  Radiation 
Therapy. 

•  Characterized  nanophosphors  suitable  for  RET,  including  Gadolinium  oxysulfide  and  Barium 
Yttrium  Fluoride. 

•  Fabricated  a  small-animal  pre-clinical  RET  imaging  box  to  enable  automated,  controlled,  RET 
experiments. 


REPORTABLE  OUTCOMES 

5  journal  publications,  including  2  first  author,  with  2  more  in  preparation;  8  conference  abstracts, 
including  3  first  author,  and  2  conference  presentations;  4  courses  taken,  including  BioE222:  Molecular 
Imaging,  Med374:  Medical  Device  Design  (taken  after  the  grant  was  accepted,  before  it  was  funded), 
the  Comprehensive  Cancer  Training  Program,  and  SIE:  The  Stanford  Institute  for  Entrepreneurs,  and  1 
undergraduate  and  1  high  school  student  mentored  (again  after  the  grant  was  accepted,  but  before  it  was 
funded).  1  patent  application. 

Oral  Presentations: 

•  Carpenter,  CM.,  et  al..  “Intraoperative  Breast  Radiotherapy  guided  by  X-ray  Luminescent  Nanoparticles,” 
Presented  at  Novel  Treatment  Delivery  Techniques,  American  Society  for  Therapeutic  Radiology  and  Oncology, 
San  Diego,  CA,  Nov.  2010. 

•  Carpenter,  CM.  “Development  of  An  X-Ray/Optical  Luminescence  Imager  for  Improved  X-Ray  Contrast 
Sensitivity”,  Presented  at  the  52nd  Annual  Meeting  of  the  AAPM,  Philadelphia,  PA,  July,  2010. 

Journal  Publications: 
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•  Carpenter  CM.,  et  al.  “Limited-Angle  X-ray  Luminescence  Tomography:  Methodology  and  Feasibility  Study” 
Phys  Med  Biol,  2011. 

•  Carpenter  CM.,  et  al.  “X-ray /Optical  Luminescence  Imaging;  Characterization  of  Experimental  Conditions,”  Med 
Phys,  37(8)  2010. 

•  Sun  C,  Pratx  G,  Carpenter  CM,  Liu  H,  Cheng  Z,  Gambhir  SS,  Xing  L.  Synthesis  and  Radioluminescence  of 
PEGylated  Eu3+-doped  Nanophosphors  as  Bioimaging  Probes  Advanced  Materials,  in  press. 

•  Sun  C,  Carpenter  CM,  Pratx  G,  Xing  L.  Facile  synthesis  of  amine  functionalized  Eu3+-doped  La(OH)3 
nanophosphors  for  bioimaging  Nanoscale  Research  Letters,  6(24)  2011. 

•  Pratx,  G.,  Carpenter,  CM.,  et  al.,  “Tomographic  Molecular  Imaging  using  X-ray  Excitable  Nanoparticles.,  Opt.  Lett., 
35(20)  3345-7,  2010. 


CONCLUSION 

The  first  funding  period  in  this  grant  has  developed  the  research  infrastructure  for  Radioluminescence 
Tomography.  This  has  resulted  in  a  fully  functioning  system  that  may  perform  the  systematic  studies  in 
phantoms  and  pre-clinical  animals  that  is  the  crux  of  Aims  3  &  4  from  the  SOW.  Towards  the  goal  of 
determining  the  feasibility  for  patient  imaging,  imaging  fiber  optic  bundles  have  been  purchased  which 
are  compatible  with  the  current  system,  and  will  provide  proof-of-concept  for  the  mammographic  and 
needle-biopsy  realizations  of  this  modality.  Future  work  for  Aims  3&4  will  involve  a  more 
sophisticated  phantom  study  to  evaluate  system  temporal  and  sensitivity  performance  with  regards  to 
imperfect  background  contrast  uptake,  and  a  systematic  study  to  determine  the  feasibility  in  pre-clinical 
and  clinical  research.  In  addition,  the  first  funding  period  for  this  grant  resulted  in  much  training  for  the 
PI,  for  molecular  imaging,  small-animal  imaging,  cancer  biology,  and  device  commercialization. 

The  research  and  training  in  the  grant  are  significant  for  the  eradication  of  breast  cancer  for  the  new 
developments  in  the  potential  for  early  detection  of  breast  cancer,  for  the  treatment  of  breast  cancer  by 
providing  feedback  to  guide  Radiation  Therapy,  for  such  breast  cancer  applications  as  Intraoperative 
Radiation  Therapy,  and  for  the  coursework  which  will  enable  this  PI  to  apply  newly  learned  skills  to 
bring  technologies  to  the  clinic  and  marketplace.  It  is  our  hope  that  this  pioneering  work  will  lead  to 
advancements  in  this  new  field  that  can  be  translated  to  the  clinic. 
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1  Hybrid  x-ray/optical  luminescence  imaging:  Characterization  of  experimental 

2  conditions 

3  C.  M.  Carpenter, a)  C.  Sun,  and  G.  Pratx 

4  Department  of  Radiation  Oncology,  School  of  Medicine,  Stanford  University,  Stanford,  California  94305 

5  R.  Rao 

6  SRI  International,  Menlo  Park,  California  94025 

7  L.  Xing 

8  Department  of  Radiation  Oncology,  School  of  Medicine,  Stanford  University,  Stanford,  California  94305 

9  (Received  21  April  2010;  revised  28  May  2010;  accepted  for  publication  6  June  2010; 

10  published  xx  xx  xxxx) 

11  Purpose:  The  feasibility  of  x-ray  luminescence  imaging  is  investigated  using  a  dual-modality 

12  imaging  system  that  merges  x-ray  and  optical  imaging.  This  modality  utilizes  x-ray  activated 

13  nanophosphors  that  luminesce  when  excited  by  ionizing  photons.  By  doping  phosphors  with  lan- 

14  thanides,  which  emit  light  in  the  visible  and  near  infrared  range,  the  luminescence  is  suitable  for 

15  biological  applications.  This  study  examines  practical  aspects  of  this  new  modality  including  phos- 

16  phor  concentration,  light  emission  linearity,  detector  damage,  and  spectral  emission  characteristics. 

17  Finally,  the  contrast  produced  by  these  phosphors  is  compared  to  that  of  x-ray  fluoroscopy. 

18  Methods:  Gadolinium  and  lanthanum  oxysulfide  phosphors  doped  with  terbium  (green)  or  eu- 

19  ropium  (red)  were  studied.  The  light  emission  was  imaged  in  a  clinical  x-ray  scanner  with  a  cooled 

20  CCD  camera  and  a  spectrophotometer;  dose  measurements  were  determined  with  a  calibrated 

21  dosimeter.  Using  these  properties,  in  addition  to  luminescence  efficiency  values  found  in  the  litera- 

22  ture  for  a  similar  phosphor,  minimum  concentration  calculations  are  performed.  Finally,  a  2.5  cm 

23  agar  phantom  with  a  1  cm  diameter  cylindrical  phosphor-filled  inclusion  (diluted  at  10  mg/ml)  is 

24  imaged  to  compare  x-ray  luminescence  contrast  with  x-ray  fluoroscopic  contrast  at  a  superficial 

25  location. 

26  Results:  Dose  to  the  CCD  camera  in  the  chosen  imaging  geometry  was  measured  at  less  than  0.02 

27  cGy/s.  Emitted  light  was  found  to  be  linear  with  dose  (R:  =  1)  and  concentration  (R2  =  1).  Emission 

28  peaks  for  clinical  x-ray  energies  are  less  than  3  nm  full  width  at  half  maximum,  as  expected  from 

29  lanthanide  dopants.  The  minimum  practical  concentration  necessary  to  detect  luminescent  phos- 

30  phors  is  dependent  on  dose;  it  is  estimated  that  subpicomolar  concentrations  are  detectable  at  the 

31  surface  of  the  tissue  with  typical  mammographic  doses,  with  the  minimum  detectable  concentration 

32  increasing  with  depth  and  decreasing  with  dose.  In  a  reflection  geometry,  x-ray  luminescence  had 

33  nearly  a  430-fold  greater  contrast  to  background  than  x-ray  fluoroscopy. 

34  Conclusions:  X-ray  luminescence  has  the  potential  to  be  a  promising  new  modality  for  enabling 

35  molecular  imaging  within  x-ray  scanners.  Although  much  work  needs  to  be  done  to  ensure  bio- 

36  compatibility  of  x-ray  exciting  phosphors,  the  benefits  of  this  modality,  highlighted  in  this  work, 

37  encourage  further  study.  ©  2010  American  Association  of  Physicists  in  Medicine. 

38  [DOI:  10.1118/1.3457332] 


39  I.  INTRODUCTION 

AQ:  40  Molecular  imaging  promises  increased  sensitivity  and  speci- 

#1  „ 

41  ficity  to  disease  compared  to  traditional  anatomical  imaging 

42  modalities.  The  information  gained  from  molecular  imaging 

43  has  the  potential  to  provide  patient-specific  selection  of 

44  therapy,  improved  prediction  of  outcomes,  and  increased 

45  treatment  efficacy.  X-ray  radiography  and  computed  tomog- 

46  raphy  (CT)  are  commonly  used  anatomical  imaging  modali- 

47  ties;  however,  although  they  provide  invaluable  information 

48  in  the  clinic,  they  have  been  largely  unsuccessful  for  molecu- 

49  lar  imaging.  This  deficiency  is  due  to  their  lack  of  sensitiv- 

50  ity  to  low  concentrations  of  contrast  agents;  x-ray  imaging  is 

51  commonly  believed  to  be  many  orders  of  magnitude  less 

52  sensitive  than  opticaF  or  radionuclide4  imaging.  This  poor 


sensitivity  arises  from  the  low  x-ray  stopping  power  of  di-  ^3 
luted  contrast  agents,  which  necessitates  high  concentrations  54 
compared  to  other  imaging  modalities.  55 

Phosphors  are  well-established  materials  used  universally  56 
in  cathode  ray  tubes  and  light-emitting  diodes  for  their  abil-  57 
ity  to  emit  light  upon  excitation  by  electrons  or  photons.  58 
Phosphors  consist  of  solid-state  crystals,  which  are  typically  59 
doped  with  transition  metals  or  lanthanide  ions.  These  mate-  60 
rials  form  a  system  optimized  to  capture  higher-energy  radia-  61 
tion  and  emit  downconverted  energy  as  optical  photons.  In  62 
the  context  of  an  x-ray  scanner,  x-ray  photons  transfer  some  63 
or  all  of  their  energy  to  electrons  in  the  solid-state  crystal  64 
through  Compton  and  photoelectric  interactions.'  These  65 
high-energy  electrons  progressively  lose  energy  through  in-  66 
teractions  with  the  atoms,  leaving  a  track  of  excited  electrons  67 
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68  behind.  When  the  energy  of  the  excited  electrons  in  the  con- 

69  duction  band  is  reduced  to  approximately  two  to  three  times 

70  the  band  gap,  they  may  migrate  into  the  luminescence  center 

71  of  the  phosphor,  the  dopant,  and  recombine  with  holes  to 

72  emit  light.6  Thus,  this  amplification  process  results  in  effec- 

73  tive  quantum  efficiencies  (photons  emitted  divided  by  the 

74  photons  absorbed),  which  can  be  much  greater  than  1.  For 

75  example,  on  the  average,  6000  photons  are  produced  for  each 

76  100  keV  x-ray  photon  absorbed  in  one  particular  gadolinium 

77  oxysulfide:terbium  phosphor.7  The  emitted  light  may  then  be 

78  imaged  by  sensitive  optical  detectors. 

79  This  paper  investigates  the  use  of  nanosized  inorganic 

80  phosphors7  as  potential  biological  contrast  agents  for  medical 

81  imaging  in  a  combined  x-ray/optical  instrument.  The  emis- 

82  sion  from  this  contrast  agent  is  evaluated  to  determine  the 

83  practicality  of  this  new  modality.  The  implications  for  this 

84  x-ray  activated  contrast  agent  are  discussed  with  regard  to  its 

85  potential  to  enable  molecular  imaging  during  fluoroscopy, 

86  x-ray  CT,  or  projection  x-ray  imaging. 


87  II.  METHODS 


Fig.  1.  (a)  Schematic  of  the  experimental  setup  including  the  CCD  camera, 
the  x-ray  source  and  detector,  and  the  sample,  (b)  Photograph  of  the  imaging 
setup. 


A  cone  beam  computed  tomography  (CBCT)  system  121 
(Acuity,  Varian  M  edical  Systems,  Palo  Alto,  CA)  was  used  to  122 
irradiate  the  sample.  This  system  allows  both  CBCT  and  123 
fluoroscopy  at  various  tube  voltages  between  70  and  125  kV.  124 
This  system  was  operated  in  fluoroscopy  mode  to  enable  125 
continuous  irradiation  from  a  constant  geometry.  126 


88  II.A.  Phosphor  fabrication 

89  Trivalent  europium  (Eu)  or  terbium  (Tb)  activated  gado- 

90  linium  oxysulfides  (GOSs)  or  lanthanum  oxysulfides  (LOSs) 

91  were  synthesized  using  appropriate  rare  earth  nitrates 

92  (99.99%  pure)  with  two  standard  methods:  The  gel-polymer 

93  combustion  process  and  the  combination  capping  process,6 

94  respectively.  After  preparation,  samples  were  heat  treated  at 

95  500-600  °C  for  1-3  h  to  aid  the  migration  of  the  dopant 

96  into  the  crystal  lattice  structures.  Next,  the  powders  were  ball 

97  milled  with  glass  beads  (10  yim)  in  the  presence  of  the  ap- 

98  propriate  surfactant  for  2-3  h. 


99  II. B.  Spectroscopy  and  imaging  of  phosphor 

100  characteristics 

101  To  facilitate  spectroscopy  and  imaging  for  the  analysis  of 

102  the  properties  of  the  phosphors,  dry  phosphor  was  placed  in 

103  plastic  test  tubes.  For  spectroscopy,  the  distal  end  of  a 

104  10  m-400  yu,m  optical  fiber  was  placed  in  contact  with  the 

105  side  of  the  test  tube,  while  the  proximal  end  was  attached  to 

106  a  spectrophotometer  (J az,  Ocean  Optics,  Dunedin,  FL), 

107  which  was  operated  from  the  x-ray  control  room.  Optical 

108  emission  spectra  across  the  visible  and  near  infrared  (NIR) 

109  range  were  acquired  with  the  spectrasuite  (Ocean  Optics) 

110  software  package.  Imaging  was  performed  with  a 

111  512x512  pixel  backi  1 1  umi  nated  CCD  camera  (CCD  tem- 

112  perature  maintained  at  -70  °C)  with  a  FI. 4  imaging  lens, 

113  exposure  times  varying  from  <1  to  60  s,  and  the  lens  aper- 

114  ture  fully  open.  During  data  acquisition,  the  imaging  camera 

115  was  shielded  with  lead  bricks  and  placed  15-20  cm  outside 

116  the  direct  field  of  radiation  to  protect  from  ionizing  radiation. 

117  The  optics  setup  was  placed  inside  a  light-tight  box  to  el i mi¬ 
ns  nate  ambient  room  light.  A  schematic  of  this  imaging  setup  is 

119  shown  in  Fig.  1(a),  and  a  photograph  of  the  setup  during 

120  experimental  imaging  is  shown  in  Fig.  1(b). 


II. C.  Concentration  evaluation  127 

Minimum  detectable  concentrations  were  calculated,  as-  128 
suming  a  superficial  location  of  the  phosphor  contrast  agent  129 
(i.e.,  no  signal  loss  due  to  tissue  absorption).  Including  light  130 
detection  losses  L,  the  light  detected  is  131 

0>  =  r  D  c-L,  (1)  132 

where  Y  is  the  emitted  light  efficiency,  D  is  the  dose,  and  cis  133 
the  concentration  of  the  phosphor.  From  Ref.  8,  the  emitted  134 
light  efficiency  (in  a  lanthanum  oxysulfide:terbium  phos-  135 
phor)  is  1.39 xlO15  optical  photons/ (Gy  mg).  We  assumed  136 
99%  detection  loss  due  to  suboptimal  optical  collection  ge-  137 
ometry.  Signal  below  a  signal-to-noise  ratio  (SNR)  of  10  138 
(assuming  shot-noise  limited  detection)  was  assumed  to  be  139 
too  low  to  detect.  To  calculate  the  molar  concentration,  we  140 
assumed  a  spherical  10  nm  diameter  nanoparticle  consisting  141 
of  hexagonal-structured  phosphors  with  lattice  constants  of  142 
a=4.046  A  and  c=6.951  A  (Ref.  9)  and  density  of  143 
5.5  g/cm3.10  144 

II. D.  Optical  phantom  fabrication  and  imaging  145 

A  small-animal  sized  tissue-simulating  phantom  was  146 
fabricated  for  this  study.  The  cylindrical  phantom  147 
measured  2.5  cm  in  diameterx4.5  cm  in  height,  with  a  148 
1x2.5  cm2  cylindrical  inclusion.  The  phantom  was  made  149 
from  1%  agar  with  homogeneous  optical  properties  using  150 
titanium  oxide  for  scatter  and  India  Ink  for  absorption  using  151 
methods  common  to  diffuse  optical  phantoms.11  The  optical  152 
properties  were  determined  by  a  previously  established  153 
system12  to  be  /j,a- 0.0025  mm-1  and  yu,s'=0.77  mm-1  at  154 
630.  Micrometer-sized  GOS:Eu  phosphor  particles  were  155 
added  to  the  inclusion  at  a  concentration  of  10  mg/ml,  and  no  156 
phosphor  was  added  to  the  background.  This  phantom  was  157 
imaged  with  an  electron  multiplied  (EM)-CCD  (Flamamatsu  158 
ImageEM  9100-13,  Flamamatsu,  Japan)  with  a  512  159 
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Fig.  2.  (a)  Radiation  dose  on  detector  for  three  different  configurations:  (i)  Direct  irradiation,  (ii)  direct  irradiation  with  lead  shielding  (5  mm  on  top,  50  mm 
on  each  side),  and  (iii)  indirect  irradiation  with  lead  shielding,  (b)  Radiation  events  per  area  per  dose. 


160  x512  pixel  sensor  cooled  at  -70  °C,  with  a  FI. 4  imaging 

161  lens,  with  exposure  of  times  less  of  2  s,  gain  at  half- 

162  maximum,  and  the  lens  aperture  fully  open. 

163  III.  RESULTS 

164  III.A.  CCD  interaction  with  ionizing  radiation 

165  To  ensure  that  the  CCD  camera  would  not  be  damaged 

166  from  the  ionizing  radiation,  an  ionization  chamber  (PTW 

167  Farmer  30010,  PTW  Freiberg  GMBH,  Germany)  was  placed 

168  in  the  radiation  field  to  determine  the  dose  to  air  in  the  vi- 

169  cinity  of  the  CCD.  Several  locations  along  the  patient  bed 

170  were  measured  to  determine  dose  rate,  as  indicated  in  Fig. 

171  2(a).  Under  direct  radiation,  the  dose  was  nearly 

172  300  yu,Gy/s;  this  dose  rate  was  reduced  by  two  orders  of 

173  magnitude  by  placing  lead  around  the  chamber,  and  further 

174  by  twofold,  to  1.5  ytGy/s,  by  moving  the  chamber  15  cm 

175  out  of  the  radiation  field.  This  rate  deposits  dose  well  below 

176  levels  that  would  likely  damage  the  CCD. 

177  To  investigate  camera  damage  further,  we  investigated  the 

178  lingering  effects  of  radiation  on  the  CCD.  The  interaction  of 

179  an  x-ray  photon  with  the  CCD  camera  appears  in  the  image 

180  as  a  bright  pixel  at  near-maximum  intensity.  These  hot  pixels 

181  appear  similar  to  cosmic  ray  interactions,  which  are  common 

182  with  CCD  cameras.  We  investigated  the  incidence  of  these 

183  events  for  a  typical  setup  using  a  CCD  to  radiation  field 

184  distance  of  15  cm  and  6  mm  of  lead  shielding  above  the 

185  camera  to  protect  from  x-ray  collimator  leakage.  Radiation 

186  events  were  determined  by  performing  an  intensity  threshold 

187  on  an  image  acquired  with  the  lens  cap  on.  It  is  clear  from 

188  Fig.  2(b)  that  the  number  of  radiation  events  is  linear  with 

189  dose,  and  thus  there  were  no  lingering  effects  from  the  ra- 

190  diation.  To  improve  image  quality,  denoising  strategies  may 

191  be  employed  utilizing  this  linearity,  such  as  an  automatic 

192  selection  of  a  hot-pixel  threshold,  which  is  dependent  on 

193  camera  dose.  Further  consideration  of  the  damage  limits  of 

194  CCD  cameras  is  given  in  Sec.  IV. 


III.B.  Phosphor  characterization  19:3 

III.B.1.  Spectral  emission  196 

A  large  body  of  knowledge  exists  on  phosphors  due  to  197 
over  a  half-century  of  study  of  optimizing  phosphors  for  198 
such  applications  as  light-emitting  diodes,  cathode  ray  tubes,  199 
and  scintillators  for  medical  imaging.  This  work  has  resulted  200 
in  a  library  of  crystals  and  dopants  from  which  one  may  201 
select  an  emission  wavelength  that  is  ideal  for  a  particular  202 
application.6  We  investigated  the  feasibility  of  GOS  phos-  203 
phors,  which  were  doped  with  either  terbium  (GOS:Tb)  or  204 
europium  (GOS:Eu),  because  of  their  absorption  K-edge  in  205 
the  diagnostic  energy  regime  at  approximately  50  keV.13  Fig-  206 
ure  3(a)  shows  the  emission  of  these  phosphors  under  100  207 
kV  x-ray  irradiation.  The  GOS:Tb  phosphor  had  a  maximum  208 
peak  emission  of  545  nm  in  green,  whereas  the  GOS:Eu  209 
phosphor  had  several  peaks  of  longer  wavelengths  in  the  210 
NIR,  including  596,  618,  627,  and  707  nm,  with  an  emission  211 
maximum  at  627  nm.  The  flexibility  enabled  by  modifying  212 
the  dopant  is  of  great  value  for  matching  the  emission  wave-  213 
length  to  a  particular  application,  such  as  the  absorption  peak  214 
of  a  phototherapeutic  drug,14  or  the  tissue  absorption  mini-  215 
mum  for  optical  imaging  in  deep  tissue.15  216 


III.B. 2.  Light  output  vs  dose  217 

To  determine  the  linearity  of  light  output  from  phosphor,  218 
GOS:Tb  was  dispersed  in  a  cuvette  containing  1%  agar  at  a  219 
concentration  of  10  mg/ml.  The  phosphor  solution  was  220 
placed  at  the  same  source-target  distance  as  the  ionization  221 
chamber.  The  x-ray  system  was  operated  in  fluoroscopy  222 
mode  and  the  tube  voltage  was  set  to  lOOkV.  The  dose  was  223 
linearly  increased  by  two  methods.  First,  the  tube  current  224 
was  increased  from  5  to  20  mA  with  a  constant  tube  voltage.  225 
Second,  the  tube  voltage  and  current  remained  constant,  and  226 
dose  was  linearly  increased  by  adjusting  the  irradiation  time.  227 
Images  were  acquired  with  a  CCD  camera,  with  the  exposure  228 
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Fig.  3.  (a)  Emission  spectra  of  GOS  phosphors  under  x-ray  excitation.  GOS:Tb,  with  an  emission  peak  at  545  nm,  is  labeled  with  asterisks  (*),  and  GOS:Eu, 
with  emission  peaks  at  596,  618,  627,  and  707  nm.  (b)  X-ray  induced  GOS:Tb  Photon  output  vs  dose  with  dose  varied  by  increasing  tube  current,  (c)  Photon 
output  vs  dose  with  dose  varied  by  increasing  irradiation  time.  For  both  methods,  photon  output  is  linear  with  dose  shown  in  terms  of  current 
(r=l,  p<0.001,  R2=l).  (d)  Photon  output  vs  concentration  for  the  GOS:Eu  phosphor.  Concentrations  were  measured  via  microcuvette  and  dispersed  in  1% 
agar.  Photon  output  is  linear  with  concentration  (r=0.99,  p=1.9xl0~7). 


229  time  optimized.  All  measurements  were  denoised  for  hot 

230  spots  and  identical  regions  of  interest  (ROIs)  were  selected 

231  for  analysis. 

aq:  232  Both  current/dose  and  time/dose  linearity  were  confirmed 
#3  233  by  running  a  linear  correlation  analysis  (correlation  coeffi- 

234  cient  of  1.0,  p<0.001).  Figure  3(b)  shows  the  linearity  in 

235  dose  using  the  first  method,  which  showed  significant  linear- 

236  ity  (correlation  coefficient  of  1.0,  p< 0.001),  while  Fig.  3(c) 

237  confirms  linearity  with  the  second  method  (correlation  coef- 

238  ficient  of  1.0,  p< 0.001). 

239  III.B.3.  Light  output  vs  concentration 

240  To  assess  the  light  output  due  to  various  phosphor  con- 

241  centrations,  dilutions  of  GOS:Eu  phosphors  from  5  yug/ml 

242  to  10  mg/ml  were  dispersed  in  1%  agar.  Cuvettes  were 

243  placed  in  a  50  kV,  30  mA  x-ray  source  and  imaged  with  a 

244  CCD  camera.  ROIs  were  selected  to  contain  similar  areas  of 

245  the  cuvettes,  and  the  signal  was  normalized  according  to  ex- 

246  posure  time.  H ot  spots  were  removed  from  the  images  prior 

247  to  analysis. 

248  The  results  shown  in  Fig.  3(d)  demonstrate  a  strong  lin- 

249  earity  (linear  correlation  coefficient  of  0.99,  p<0.001)  with 

250  a  slope  of  1.09.  The  slight  departure  of  the  slope  from  unity 

251  is  most  likely  due  to  errors  in  selecting  identical  ROIsforthe 


different  cuvettes  and  may  have  resulted  in  the  inclusion  of  252 
the  cuvette  wall  in  the  ROI,  which  exhibited  some  light  pip-  253 
ing.  254 

III.C.  Minimum  detectable  concentrations  255 

We  calculated  the  minimum  detectable  concentration  ac-  256 
cording  to  the  methodology  outlined  in  Sec.  II  C,  for  doses  257 
ranging  from  1  cGy  (less  than  the  typical  mammographic  258 
dose)  to  20  Gy  (a  typical  dose  delivered  in  single-dose  intra-  259 
operative  radiation  therapy).  In  addition,  we  calculated  the  260 
minimum  detectable  concentration  for  several  phosphor  effi-  261 
ciencies,  scaled  according  to  that  reported  by  Kandarakis  et  262 
a/.8  (i.e.,  100%  is  equivalent  to  the  efficiency  reported).  Fig-  263 
ure  4  shows  the  minimum  detectable  concentrations  for  these  264 
scenarios.  According  to  these  calculations,  picomolar  (ng/ml)  265 
concentrations  are  detectable  (SNR  of  10)  with  a  LOS  phos-  266 
phor  for  mammographiclike  dose,  while  therapy  doses  allow  267 
femtomolar  (pg/ml)  concentrations  to  be  detected.  268 

III.D.  Contrast  comparison  between  x-ray/optical  and  269 
fluoroscopy  270 

The  recovered  contrast  between  an  inclusion  with  phos-  271 
phor  and  a  background  without  phosphor  was  investigated  to  272 


Medical  Physics,  Vol.  37,  No.  8,  August  2010 


5 


Carpenter  eta/.:  X-ray  luminescence  imaging 


5 


a  b 


Fig.  4.  M  inimum  estimated  phosphor  concentrations  calculated  from  the  literature  and  our  own  data.  Lines  represent  phosphor  luminescence  efficiencies;  the 
100%  line  is  based  on  efficiency  data  from  a  La202S  :Tb  phosphor  reported  by  K  andaraki s  etal.  (a)  M  inimum  weight  concentration  (in  mg/ml)  vs  dose  (Gy), 
(b)  M  inimum  molar  concentration  (in  M  )  vs  dose  (Gy). 


273  compare  the  contrast  differences  between  x-ray  fluoroscopy 

274  and  x-ray/optical  luminescence  imaging  in  the  small-animal 

275  imaging  phantom  described  in  Sec.  II  D.  The  phantom  was 

276  imaged  during  fluoroscopy  operation  with  the  tube  voltage  at 

277  100  kV  and  the  tube  current  at  10  mA .  The  x-ray  source  was 

278  placed  above  the  phantom  to  evenly  irradiate  the  volume. 

279  The  CCD  camera  was  placed  within  15  cm  of  the  phantom 

280  and  oriented  orthogonal  to  the  direction  of  irradiation.  The 

281  fluoroscopy  image  was  taken  simultaneously.  It  is  important 

282  to  note  that  gadolinium  is  a  strong  absorber  of  100  keV  x-ray 

283  energy,  with  a  mass  attenuation  coefficient  of  3.109  cm2/g 

284  (compared  to  common  x-ray  contrast  agents  such  as  barium 

285  at  2.196  cm2/g  and  iodine  at  1.942  cm2/g).  Since  gado- 

286  I i ni um  has  a  higher  mass  attenuation  coefficient  for  x-ray 

287  photons  than  water  (mass  attenuation  coefficient  of 

288  0.1707  cm2/g),  it  should  exhibit  slight  contrast. 

289  The  images  from  the  phantom  are  shown  in  Fig.  5.  Figure 

290  5(a)  shows  a  white  light  image  taken  by  the  CCD  camera 

291  with  background  illumination  from  the  room  lights.  The  cor- 

292  responding  fluoroscopy  image  is  shown  in  Fig.  5(b).  The 

293  phosphor  inclusion  is  indicated  by  the  red  arrow  and  shows 

294  slight  increased  x-ray  absorption  (the  smaller  dark  circle  is  a 

295  bolt  hole  in  the  aluminum  optical  table).  Figure  5(c)  shows  a 

296  raw  optical  image  taken  of  the  x-ray  luminescent  phosphor. 

297  This  image  is  overlaid  on  the  white  light  image  in  Fig.  5(d). 

298  The  contrast  between  the  inclusion  and  the  background  is 

299  very  slight,  at  0.6%  for  the  fluoroscopy  image,  while  it  is 

300  over  260%  for  the  luminescence  image,  in  addition,  the 

301  signal-to-noise  ratio  for  the  phosphor  emission  was  23  vs  2.4 

302  for  the  fluoroscopy  image. 

303  IV.  DISCUSSION 

304  We  found  that  dose  distributed  to  the  shielded  camera  was 

305  measured  at  less  than  3  yctGy/s  when  the  camera  was  posi- 

306  tioned  at  the  isocenter  and  the  x-ray  tube  voltage  was  100 

307  keV  with  the  tube  current  at  20  mA.  Although  this  dose  is 

308  low,  estimating  a  damage  threshold  is  difficult  for  CCDs  be- 

309  cause  damage  is  design/manufacturer-specific  and  is  depen- 

310  dent  on  environmental  conditions  (for  a  more  thorough  over- 

311  view,  see  Ref.  16).  It  is  well  recognized  that  the  largest 


radiation  threat  to  the  operation  of  a  CCD  is  the  bombard-  312 
ment  by  highly  energetic  heavy  particles,  such  as  protons  and  313 
neutrons.  These  particles  contribute  most  to  CCD  damage  314 
through  impact  displacements  of  silicon  atoms  which  create  315 
semipermanent  energy  traps.  These  traps  create  energy  levels  316 
which  can  increase  J  ohnson  noise  (via  promoting  valence  317 
band  electrons  to  the  conduction  band),  create  spurious  noise  318 
when  trapped  electrons  are  released,  and  alter  the  operation  319 
of  transistor  gates  by  altering  their  flat-band  voltage  (for  a  320 
more  thorough  review,  see  Refs.  17  and  18).  Although  the  321 
probability  of  creating  protons  and  neutrons  is  extremely  low  322 
at  diagnostic  x-ray  energies  studied  here,  it  is  relevant  for  323 
therapeutic  energies  in  the  M  V  range.  Cumulative  doses  are  324 
also  important  because  of  the  increased  probability  for  a  325 
high-energy  photon  interaction,  it  has  been  reported  that  total  326 


Fig.  5.  Agar  imaging  phantom  with  embedded  phosphors  and  tissue  emu¬ 
lating  optical  properties,  (a)  White-light  optical  image,  (b)  Projection  fluo¬ 
roscopy  image  (note  the  distinction  between  the  phosphor  inclusion- 
indicated  by  the  red  arrow—  around  300  units,  compared  to  the  black  circle 
caused  by  a  screw  hole  in  the  optical  table  supporting  the  phantom),  (c) 
Optical  emission  from  the  phantom,  (d)  Overlay  of  the  white  light  image  (a) 
and  the  light  emission  (c). 
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327  doses  above  about  25  and  100  Gy  are  considered  thresholds 

328  for  increased  noise  and  permanent  damage,  respectively19 

329  (though  again,  these  numbers  are  highly  camera-dependent). 

330  Thus,  it  is  important  to  keep  radiation  dose  on  the  CCD  as 

331  low  as  possible  to  minimize  cumulative  dose  and  increase 

332  the  life  of  the  camera.  The  dose  deposited  in  this  experiment 

333  of  1-3  yuGy/s  is  minimal,  and  damage  due  to  radiation 

334  should  be  insignificant  in  the  lifetime  of  the  camera. 

335  The  lanthanide  dopants  examined  in  this  study  are  known 

336  to  have  extremely  sharp  peaks;  for  example,  we  measured 

337  the  full  width  at  half  maximum  (FWHM)  of  the  GOS:Tb 

338  sample  at  3  nm,  while  the  GOS:Eu  sample  had  a  FWHM  of 

339  2  nm  for  the  largest  peak  [both  shown  in  Fig.  3(a)].  This  is  in 

340  contrast  to  common  fluorophores  like  indocyanine  green, 

341  which  have  often  have  emission  peak  FWHM  of  about  60 

342  nm.  This  enables  the  possibility  of  multiplexing  contrast 

343  agents  with  very  little  cross-talk,  allowing  the  simultaneous 

344  measurement  of  several  markers  for  disease. 

345  A  major  concern  of  these  x-ray  excitable  phosphor  con- 

346  trast  agents  is  that  they  require  ionizing  dose  to  activate. 
347 Thus,  lower  concentrations  of  phosphors  will  necessitate 

348  higher  doses.  We  analyzed  the  practical  detection  limit  of 

349  phosphors  using  the  knowledge  of  concentration  and  dose 

350  linearity  and  luminescence  efficiencies  found  in  the  literature 

351  for  similar  nanophosphors.  For  detection  at  the  surface  of 

352  tissue,  x-ray  mammographic  dose  should  be  sufficient  to  al- 

353  low  the  detection  of  picomolar  (ng/ml)  concentrations  of 

354  phosphor.  This  finding  is  corroborated  with  our  experimental 

355  results  shown  in  Fig.  3.  If  the  data  in  Fig.  3(d)  are  extrapo- 

356  lated  to  10  counts/pixel,  adequate  signal-to-noise  ratio  is 

357  achieved  using  an  EM  -CCD  camera  (with  dark  noise  of  0.05 

358  counts/pixel  for  this  acquisition).  Placing  the  camera  within  a 

359  few  cm  (compared  to  20  cm  in  this  experiment)  would  result 

360  in  the  detection  of  approximately  two  orders  of  magnitude 

361  more  light.  Further,  if  the  dose  is  increased  by  an  order  of 

362  magnitude  [from  the  4  mGy  dose  administered  for  Fig.  3(d)] 

363  to  mammography  levels,  and  EM  gain  is  increased  to  the 

364  maximum,  concentrations  of  ng/ml  (picomolar)  may  be  real- 

365  ized. 

366  At  deeper  depths,  however,  light  attenuates  rapidly;  for 

367  example,  in  breast  tissue,  a  depth  of  3  cm  would  attenuate 

368  detected  light  by  approximately  three  orders  of  magnitude.20 

369  This  would  effectively  decrease  phosphor  sensitivity  during 

370  mammographic  dose  to  nanomolar  concentrations.  In  a  ra- 

371  diation  oncology  setting,  however,  high  doses  are  delivered 

372  to  treat  disease  (such  as  the  20  Gy  single-fraction  therapy, 

373  which  is  used  during  intraoperative  radiation  therapy21).  This 

374  technique  could  aid  a  surgeon  and  radiation  oncologist  to 

375  identify  the  distribution  of  disease  around  a  tumor  margin, 

376  such  as  during  breast  lumpectomy.  In  this  case,  the  patient 

377  would  be  injected  with  a  phosphor  contrast  agent  before  or 

378  during  surgery,  depending  on  the  pharmacokinetics  of  the 

379  tracer.  The  tracer  could  be  imaged  during  the  first  1%  or  10% 

380  (e.g.,  0.2  or  2  Gy)  of  the  radiotherapy  treatment  dose,  which 

381  would  provide  the  clinicians  with  more  confidence  about  the 

382  treatment  dose  or  volume,  or  enable  an  adjustment  to  the 

383  dose  distribution. 
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The  other  area  of  concern  with  these  particles  is  biologi-  384 
cal  compatibility.  This  is  an  issue  for  all  nanoparticle  sys-  385 
terns,  and  much  work  is  being  done  to  develop  strategies  to  386 
ensure  stability  and  compatibility.22  In  fact,  multiple  groups  387 
have  successfully  used  upconversion  phosphors  in  biological  388 
small-animal  models.23,24  The  increased  interest  in  phosphors  389 
should  aid  in  the  rapid  advancement  in  biocompatibility,  390 
which  will  aid  this  modality.  391 

We  found  that  the  contrast  to  background  ratio  was  over  392 
2.5  orders  of  magnitude  higher  for  optical  detection  of  the  393 
luminescent  phosphors  compared  to  x-ray  fluoroscopy.  The  394 
actual  contrast  amplification  is  much  higher  since  the  optical  395 
photons  emitted  from  the  inclusion  exhibit  extensive  scatter  396 
in  the  background  and  subsequently  are  detected  by  the  397 
CCD.  These  scatter  effects  would  be  greatly  minimized  via  398 
modeling  of  the  light  propagation.  In  comparison,  x-ray  pho-  399 
ton  scatter  is  relatively  nonexistent  so  that  the  signal  contri-  400 
bution  from  the  background  originating  from  the  inclusion  is  401 
negligible.  Considering  the  photon  scatter  should  greatly  im-  402 
prove  the  contrast  of  these  phosphors  for  optical  detection  403 
compared  to  fluoroscopy.  404 

Although  this  work  demonstrated  the  potential  of  x-ray  405 
luminescence  imaging  for  imaging  a  superficial  object,  im-  406 
aging  of  lesions  centimeter  deep  should  be  possible,  with  407 
contrast- resolution  limitations  depending  on  tissues'  proper-  408 
ties,  concentration,  and  nonspecific  uptake.  The  development  409 
of  deep-tissue  x-ray  luminescence  imaging  will  require  the  410 
incorporation  of  optical  tomographic  models.  With  x-ray  lu-  411 
minescent  imaging,  the  x-ray  source  must  be  modeled  in  412 
tissue  to  give  an  accurate  description  of  dose.  There  are  413 
many  sophisticated  tools  to  model  dose,  such  as  M  onte  Carlo  414 
or  analytical  models,  which  have  been  shown  to  be  accurate  415 
[e.g.,  within  4%  (Ref.  25)]  in  biological  tissues.  Concurrent  416 
x-ray  structural  imaging  will  further  improve  these  calcula-  417 
tions.  After  dose  distribution  is  calculated,  tomographic  im-  418 
aging  may  be  performed  with  a  reconstruction  model  that  419 
uses  a  model  of  the  light  propagation  in  tissue  to  minimize  420 
the  difference  between  calculated  and  optical  measurements.  421 
This  is  very  similar  to  the  fluorescence  molecular  imaging  422 
problem.26  Once  again,  the  knowledge  of  anatomical  infor-  423 
mation  will  aid  the  optical  reconstruction  problem  by  provid-  424 
ing  structural  detail  which  may  be  used  to  improve  optical  425 
modeling27  and  reconstruction.28  426 

The  joint  use  of  x-ray  activated  phosphors  for  molecular  427 
imaging  offers  several  advantages  to  x-ray  imaging  and  to  428 
all-optical  fluorescent  imaging.  For  x-ray  imaging,  contrast  429 
agent  imaging  is  currently  limited  to  high  concentrations  of  430 
nonspecific  iodine  or  barium  sulfate.  Optimal  concentrations  431 
for  these  contrasts  have  been  reported  around  300-500  432 
mg/ml.29  These  high  concentrations  are  impractical  for  im-  433 
aging  biological  targets.2  The  ability  to  image  cellular  targets  434 
would  be  a  great  benefit  to  x-ray  imaging,  which,  despite  435 
being  the  most  prominent  modality  in  use  in  the  clinic  today,  436 
is  generally  limited  to  imaging  structural  anatomy.  The  use  437 
of  phosphors  combined  with  the  sensitivity  of  optical  imag-  438 
ing  allows  lower,  more  biologically  feasible  concentrations  439 
of  contrast  agents  than  is  currently  available  with  x-ray  im-  440 
aging  alone.  441 
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442  The  use  of  x-ray  activated  phosphors  offers  three  unique 

443  advantages  to  all-optical  approaches.  First,  this  dual- 

444  modality  instrument  offers  inherent  spatial  co registration  be- 

445  tween  anatomical  features  and  optical  contrasts.  This  regis- 

446  tration  is  critical  for  imaging  functional  pathology  in  deep 

447  tissue  [hence,  the  need  for  positron  emission  tomography 

448  (PET)/CT  imaging  systems].30  Next,  the  use  of  x-ray  excita- 

449  tion  eliminates  the  optical  autofluorescence  issue  in  optical 

450  imaging.  Since  the  x-ray  excitation  spectrum  is  undetectable 

451  with  photo-optical  detectors,  autofluorescence  is  avoided, 

452  which  potentially  reduces  the  detection  limit  for  low  concen- 

453  trations.  Finally,  this  technique  is  also  expected  to  have  in- 

454  creased  depth  performance  over  optical  imaging,  because  of 

455  the  high  penetration  of  x-ray  photons  in  tissues.  X-ray  pho- 

456  tons  have  nearly  two  orders  of  magnitude  lower  effective 

457  attenuation  coefficient  compared  to  optical  photons;  this  op- 

458  portunity  offers  the  potential  to  use  clinically  available  in- 

459  strumentation  as  an  external  source. 


460  V.  CONCLUSIONS 

461  This  study  focused  on  the  instrumentation  and  material 

462  feasibility  of  inorganic  downconversion  phosphors  toward 

463  the  realization  of  x-ray  molecular  imaging.  Significant  recent 

464  advances  in  PET,4  optical  imaging,11  magnetic  resonance 

465  imaging,12  and  to  a  lesser  extent,  single  positron  emission 

466  computed  tomography,  and  ultrasound  have  invigorated  the 

467  search  for  disease- specific  protein  receptors  that  may  be  tar- 

468  geted  with  imaging  agents.  This  approach  has  already  been 

469  applied  to  numerous  pathologies  to  identify  atherosclerosis 

470  and  thrombosis,31  to  determine  treatment  efficacy  via  apop- 

471  tosis  markers,32  to  identify  cancer,  and  to  monitor  cellular 

472  activity.  The  incorporation  of  these  markers  into  x-ray  imag- 

473  ing  may  have  significant  impact  on  medical  imaging. 

474  In  this  work,  we  demonstrate,  for  the  first  time  to  our 

475  knowledge,  the  feasibility  of  using  inorganic  phosphors  to 

476  enable  optical  detection  under  x-ray  irradiation,  which  may 

477  enable  x-ray  molecular  imaging.  We  first  investigated  the 

478  practical  feasibility  of  operating  a  CCD  within  an  x-ray  ex- 

479  citation  field  at  clinically  relevant  energies,  taking  into  con- 

480  sideration  noise  and  potential  damage.  We  found  that  the 

481  dose  distribution  to  air  was  sufficiently  low  to  prevent  dam- 

482  age  during  operation.  Additionally,  the  noise  on  the  CCD  due 

483  to  incoming  high-energy  irradiation  is  manageable.  We  dem- 

484  onstrated  the  ability  of  several  phosphors  to  emit  light  in  the 

485  optical  regime  under  x-ray  excitation.  These  phosphors 

486  should  be  effective  for  tuning  light  output  to  a  specific  ap- 

487  plication.  We  found  that  the  light  output  was  linearly  propor- 

488  tional  to  both  dose  and  concentration.  Future  work  will  focus 

489  on  quantitative  imaging.  We  calculated  minimum  detectable 

490  concentrations  based  on  these  data  and  values  found  in  lit- 

491  erature;  these  concentrations  are  sufficient  for  certain  bio- 

492  logical  imaging  applications.  Finally,  we  demonstrated  the 

493  potential  of  inorganic  phosphors  to  image  lower  concentra- 

494  tions  than  is  possible  with  x  ray  alone.  We  found  a  430  times 

495  improvement  in  contrast  recovery  for  optical  detection  corn- 

496  pared  to  fluoroscopic  detection.  This  improvement  is  ex- 

497  pected  to  be  greater  with  modeling  of  photon  propagation 
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and  imaging.  We  envision  hybrid  x-ray/optical  imaging  may  498 
have  significant  application  in  the  detection  and  diagnosis  of  499 
disease,  especially  during  image-guided  intervention.  500 
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Abstract 

X-ray  luminescence  tomography  (XLT)  has  recently  been  proposed  as  a  new 
imaging  modality  for  biological  imaging  applications.  This  modality  utilizes 
phosphor  nanoparticles  which  luminesce  near-infrared  light  when  excited  by 
x-ray  photons.  The  advantages  of  this  modality  are  that  it  uniquely  combines 
high  sensitivity  of  radioluminescent  nanoparticles  and  high  spatial  localization 
of  collimated  x-ray  beams.  Currently,  XLT  has  been  demonstrated  using 
x-ray  spatial  encoding  to  resolve  the  imaging  volume.  However,  there 
are  applications  where  the  x-ray  excitation  may  be  limited  by  geometry, 
where  increased  temporal  resolution  is  desired,  or  where  a  lower  dose  is 
mandatory.  This  paper  extends  the  utility  of  X  LT  to  meet  these  requirements 
by  incorporating  a  photon  propagation  model  into  the  reconstruction  algorithm 
to  recover  dimensions  remaining  in  an  x-ray  limited-angle  (LA)  geometry. 
This  enables  such  applications  as  image-guided  surgery,  where  the  ability  to 
resolve  lesions  at  depths  of  several  centimeters  can  be  the  key  to  successful 
resection.  The  hybrid  x-ray/diffuse  optical  model  is  first  formulated  and 
then  demonstrated  in  a  breast-sized  phantom,  simulating  a  breast  lumpectomy 
geometry.  Both  numerical  and  experimental  phantoms  are  tested,  with  lesion- 
simulating  objects  of  various  sizes  and  depths.  Results  show  localization 
accuracy  with  median  error  of  2.2  mm,  or  4%  of  object  depth,  for  small 
2-14  mm  diameter  lesions  positioned  from  1  to  4.5  cm  in  depth.  This  compares 
favorably  with  fluorescence  optical  imaging,  which  is  not  able  to  resolve  such 
small  objects  at  this  depth.  The  recovered  lesion  size  has  lower  size  bias  in 
the  x-ray  excitation  direction  than  the  optical  direction,  which  is  expected  due 
to  the  increased  optical  scatter.  However,  the  technique  is  shown  to  be  quite 
invariant  in  recovered  size  with  respect  to  depth,  as  the  standard  deviation  is  less 
than  2.5  mm.  Sensitivity  is  a  function  of  dose;  radiological  doses  are  found  to 
provide  sufficient  recovery  for  /zg  mb1  concentrations,  while  therapy  dosages 
provide  recovery  for  ng  mh1  concentrations.  Experimental  phantom  results 
agree  closely  with  the  numerical  results,  with  positional  errors  recovered  within 
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8.6%  of  the  effective  depth  for  a  5  mm  object,  and  within  5.2%  of  the  depth 
for  a  10  mm  object.  Object-size  median  error  is  within  2.3%  and  2%  for  the 
5  and  10  mm  objects,  respectively.  For  shallow-to-medium  depth  applications 
where  optical  and  radio-emission  imaging  modalities  are  not  ideal,  such  as  in 
intra-operative  procedures,  LAXLT  may  be  a  useful  tool  to  detect  molecular 
signatures  of  disease. 

(Some  figures  in  this  article  are  in  colour  only  in  the  electronic  version) 


L  Introduction 

Imaging  plays  a  vital  role  in  the  management  of  cancer  care,  for  detection,  staging,  intervention, 
and  monitoring  of  treatment  response.  Despite  its  ubiquitous  use  elsewhere,  the  role  of  imaging 
in  surgery  is  limited,  as  it  is  dominated  by  C-arm  fluoroscopy  and  optical  endoscopy.  These 
tools  are  appropriate  for  visualizing  tissue  structure,  yet  are  limited  in  their  sensitivity  to 
microscopic  disease.  This  limitation  affects  such  procedures  as  surgical  breast  lumpectomy, 
as  many  studies  have  found  that  surgeons  are  unable  to  remove  all  tumor  tissue  present  in  the 
surgical  field  (for  example,  Gibson  etal  (2001)  identified  residual  tumor  in  55%  of  the  cases. 
The  risks  of  local  failure  are  high,  as  local  failure  often  leads  to  distant  metastasis  (Fortin  etal 
1999).  Thus,  there  is  a  need  for  tools  to  provide  surgeons  with  more  sensitive,  more  specific 
image  guidance. 

This  need  may  be  fulfilled  with  molecular  imaging,  which  promises  to  image  molecular 
and  cellular  processes,  and  may  allow  the  early  identification  of  disease  or  status  of  disease 
progression  and  treatment  (Weissleder  and  Pittet  2008).  Developing  these  tools  for  the 
operating  room  would  aid  a  physician  during  an  intervention,  by  allowing  the  clinician  to 
identify  near-microscopic  regions  of  disease,  such  as  at  the  tumor  margin.  Ideally,  this  tool 
would  be  able  to  image  at  a  depth  of  several  centimeters,  so  that  disease  buried  beneath 
the  superficial  layers  could  be  identified.  Several  potential  applications  for  this  technology 
could  be  in  removing  occult  disease  in  breast  (Tanaka  etal  2006,  Alex  and  Krag  1993),  brain 
(Stummer  etal  2008)  and  hepatic  tumors  (Torzilli  etal  1999),  where  imaging  is  currently 
being  incorporated  into  the  clinic,  and  new  innovations  may  be  readily  translated. 

This  paper  develops  and  demonstrates  a  novel  x-ray  luminescence  tomographic  (XLT) 
method  that  is  uniquely  suited  for  image-guided  surgical  applications.  This  method,  Limited- 
AngleX  LT  (LAX  LT ),  uti lizes  a  photon  propagation  model  to  enableX  LT  forsurgical  guidance, 
where  X  LT's  advantages  are  the  clearest  for  translation  into  the  clinic.  X  LT  has  been  recently 
introduced  (Carpenter  etal  2010)  and  demonstrated  in  simulation  and  in  phantoms  (Pratx  etal 
2010a,  2010b).  This  imaging  modality  utilizes  nano-sized  phosphors  which  emit  optical  near- 
infrared  light  upon  x-ray  excitation  (Chen  2008,  Sun  etal  2010).  Attaching  these  phosphors 
to  molecular  probes  (e.g.  antibodies  and  peptides)  that  target  molecular  markers  specific  to 
tumors,  such  as  angiogenesis  markers  like  epidermal  growth  factor  receptor  (Sokolov  etal 
2003),  or  a„^3-integrin  (H aubner  eta/  2001,  Chen  etal  2004)  expression,  could  allow  the 
surgeon  to  differentiate  between  normal  and  cancerous  tissue.  X  LT  has  several  advantages  to 
current  molecular-sensitive  imaging  modalities:  emission  imaging  techniques,  such  as  gamma 
cameras,  are  limited  in  their  ability  to  discriminate  depth  due  to  the  limited  angles  that  may  be 
imaged  during  surgical  procedures  (Barrett  1990);  optical  imaging,  on  the  other  hand,  has  the 
ability  to  provide  depth  localization,  and  is  currently  under  investigation  for  surgical  guidance 
(Tanaka  etal  2006,  Stummer  etal  2008,  Roberts  etal  2010),  yet  is  limited  in  its  ability  to 
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image  deeper  than  ~1  cm  (Kepshire  et  al  2007).  Depth  is  important  to  discern  to  determine 
occult  lesions  lying  under  the  superficial  layer,  and  to  determine  the  feasibility  of  surgical 
removal  of  a  lesion. 

XLT  utilizes  the  extremely  low  scatter  of  x-rays  compared  to  optical  fluorescence  imaging 
to  enable  higher  spatial  resolution.  A  thin  pencil-beam  of  collimated  x-rays  may  be  maintained 
while  the  x-rays  propagate  through  tissue  of  several  cm;  this  spatial  localization  is  in  contrast 
to  optical  excitation,  which  is  highly  attenuated  and  scattered  (O'Leary  eta/  1995).  By  rotating 
the  x-ray  (or  similarly  the  phantom)  to  cover  all  angular  projections,  the  resolving  power  is 
limited  merely  by  the  width  of  the  beam  (up  to  the  diffraction  of  the  x-ray).  A  numerical 
analysis  demonstrated  that  2.25  cm  deep  objects  as  small  as  1  mm  (using  a  1  mm  beam  width) 
with  a  nanoparticle  concentration  of  0.4  pM  could  be  resolved;  increasing  dose  increased  the 
sensitivity  (Pratx  et  al  2010a).  However,  there  are  applications  where  the  x-ray  excitation 
may  be  limited  by  geometry,  where  increased  temporal  resolution  is  desired,  or  where  a  lower 
dose  is  mandatory;  one  such  application  is  intraoperative  breast  cancer  lumpectomy,  where 
it  may  not  be  possible  or  desirable  to  irradiate  over  the  full  projection  space.  In  these  cases, 
it  would  be  beneficial  to  irradiate  over  a  limited  projection  space,  and  use  the  ability  of  the 
optical  detectors  to  resolvethe  remaining  dimensions.  Thistechniquecould  also  have  utility  in 
decreasing  dose  to  the  tissue,  as  fewer  irradiation  beamlets  are  needed  to  resolvethe  volume. 

This  paper  develops  a  reconstruction  methodology  for  utilizing  XLT  to  perform  depth- 
resolved  imaging  in  a  geometry  appropriate  for  tumor-resection  applications.  This  method 
develops  a  hybrid  x-ray/optical  reconstruction,  which  allows  XLT  spatial  encoding  in 
a  limited-angle  geometry,  and  diffuse  optical  spatial  discrimination  for  the  remaining 
dimensions.  This  technique  augments  that  of  Pratx  eta/  (2010a),  who  encoded  all  spatial 
dimensions;  such  a  technique  is  more  suitable  for  such  applications  as  small-animal  imaging. 
The  advantage  of  this  new  approach  is  that  enables  X  LT  in  surgical  applications  such  as  breast 
or  brain  excision,  and  may  reduce  dose.  The  performance  of  this  technique  is  examined  in 
both  numerical  and  experimental  phantoms  for  various  object  sizes  and  positions,  within  a 
geometry  that  mimics  breast  and  brain  intraoperative  geometries. 


Z  Methods 

2.1.  Experimental  setup 

Theequipment  used  for  this  study  consisted  of  an  x-ray  radiation  sourceto  excite  the  phosphors, 
and  an  optical  detector  to  sample  the  photon  fluence.  During  acquisition,  the  radiation  source 
is  collimated  into  a  thin  slice  as  described  by  Pratx  et  al  (2010a)  to  excite  a  plane  shaped 
volume.  An  optical  camera  samples  the  emitted  light.  This  schematic  is  shown  in  figure  1(a). 
The  experimental  setup  is  shown  in  figure  1(b). 

The  x-ray  source  used  for  the  measurements  in  this  paper  was  a  50  kVp  x-ray  superficial 
unit  (Pantak  Therapax-150,  Elimpex,  AT)  with  a  10  cm  exit-diameter  cone  applicator.  This 
cone  was  placed  17.5  cm  away  from  surface  of  the  phantom.  T  he  beam  was  coll  imated  to  1  mm 
wide  (verified  optically)  by  carefully  positioned  50  mm  thick  lead  bricks.  A  high-sensitivity 
EM -CCD  camera  (Pro-EM,  Princeton  Instruments,  NJ)  with  an  F/1.4  lens  was  positioned 
~20  cm  away  from  the  surface  of  the  phantom.  This  distance  was  chosen  to  minimize  x- 
ray  photon  noise  (Carpenter  eta/  2010).  As  an  alternative,  optically  clear  leaded  acrylic  or 
leaded  glass  could  be  placed  between  the  camera  and  the  sample  and  used  to  reduce  x-ray 
noise  on  the  CCD  and  allow  the  camera  to  be  placed  closer  to  the  sample  to  collect  more 
light. 
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Figure  1  (a)  Schematic  of  LAXLT  while  imaging  a  phantom  with  a  single  spherical  object, 
(b)  Experimental  setup  used  in  this  study.  A  black-box  was  used  to  enclose  the  phantom  (not 
shown)  and  eliminate  ambient  light  from  the  experiment. 


2.2.  Image  formation 

2.2.1.  X-ray  luminescence  emission  forward  model.  To  determine  the  concentration  of 
nanophosphors,  the  radiation  (x-ray  and  optical)  must  be  modeled.  The  emission  of  x-ray 
excitable  nanophosphors  is  linearly  dependent  on  the  dose  imparted  to  the  tissue  (Carpenter 
et  al  2010),  d  (units  of  Gy),  the  fractional  efficiency  of  the  phosphor  in  converting  ionizing 
energy  to  optical  emission,  r,  and  the  concentration,  c  (units  of  mg  mh1).  The  luminescent 
photon  density  from  the  nanophosphors,  $,  due  to  an  incident  radiation  beam  is  therefore 

<£>  =  Fdc.  (1) 

Determining  the  total  ionization  energy  imparted  to  the  tissue  (dose)  is  a  procedure  that 
requires  calibration  to  incorporate  the  properties  of  the  radiation  emitted  by  thex-ray  system. 
This  calibration  is  system  specific,  taking  into  account  tube  potential,  geometry,  x-ray  tube 
target,  and  filter  material;  these  factors  taken  together  form  an  x-ray  spectrum,  known  as 
the  beam  quality.  The  dose  at  depth  is  determined  using  measurements  from  a  calibrated 
ionization  chamber  in  a  phantom  and  composed  into  a  look-up  table,  the  percent  depth  dose 
(PDD)  curves.  This  system- wide  calibration  is  performed  periodically  (M  a  etal  2001).  Using 
the  PDD  curves,  dose  at  a  specific  depth  in  the  tissue  can  be  determined  by  knowing  the 
source-to-surface  distance  between  the  x-ray  tube  and  the  tissue.  This  method  can  have 
high  quantitative  accuracy  of  1-2%  (M  unck  af  Rosenschold  etal  2008).  Another  method  to 
accurately  determinedoseisthrough  M  onteCarlo  methods,  which  model  the  system,  including 
the  above  factors  and  also  including  patient  anatomy.  This  method  can  calculate  dose  with 
high  accuracy  as  long  as  comprehensive  modeling  of  beam  quality  is  performed  (Verhaegen 
etal  1999).  In  this  study,  we  used  the  PDD  curves  to  determine  dose. 


2.2.2.  Diffuse  optical  forward  model.  Images  acquired  at  the  tissue  surface  are  input  into 
a  photon  propagation  model  to  determine  the  phosphor  distribution.  The  images  from  the 
CCD  camera  are  first  processed  to  remove  x-ray  noise  using  a  simple  gradient-threshold 
algorithm,  and  then  input  into  the  algorithm  as  the  data.  The  propagation  of  optical  light  can 
be  approximated  by  the  lossy  photon  diffusion  equation  (DE)  (Arridge  eta/  1993),  which 
yields  the  photon  density  in  tissue.  The  DE  is  valid  for  many  soft  human  tissues,  including 
the  breast,  lung,  prostate,  brain,  etc  (Cheong  etal  1990).  Following  excitation  from  x-ray 
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radiation,  the  time- independent  luminescence  photon  density  emitted  from  the  nanophosphors 
is 


O(r)  =  -V-D(r)V«r)  +  /ia(r)0(r)  (2) 

where  4>[r)  is  the  photon  fluence  at  position  r,  in  units  of  photons  per  area  per  time,  and  4>(r) 
is  the  photon  density,  in  units  of  photons  per  volume  per  time.  Photon  propagation  is  affected 
by  the  absorption  and  diffusion  coefficients  of  the  tissue,  /j.a  and  D,  respectively,  which  are 
dependent  on  wavelength.  Thediffusion  coefficient,  D,  isdefined  as  D  =  3  1  ,  where  n's  is 

the  reduced  scattering  coefficient  of  the  emitted  photons.  A  type  I II  boundary  condition  is  used 
to  model  the  photon  fluence  at  the  boundary,  -DV0  •  ft  =  a<t>,  where  a  defines  the  internal 
reflection  of  the  light  at  the  tissue  boundary  due  to  the  index  of  refraction  mismatch  between 
tissue  and  air  (Schweiger  eta/  1995,  Aronson  1995),  and  the  unit  vector  n  is  normal  to  the 
surface  of  the  phantom.  Because  no  unique  solution  exists  for  (2)  with  arbitrary  boundaries, 
equation  (2)  is  approximated  with  the  finite  element  method  (FEM  )  (A  rridge  eta/ 1993).  This 
problem  is  similar  to  the  diffuse  optical  fluorescence  model  introduced  by  J  iang  (1998),  and 
is  adapted  here. 

As  described  by  J  iang,  the  photon  emission  may  be  approximated  with  the  FEM  by 

[A]{0>  =  {*}  (3) 

where  A  i s  the  FEM  approximation  of  the  physics  of  photon  propagation  (the  right-hand  sideof 
equation  (2))  and  bis  the  approximation  to  the  light  source  (the  left-hand  sideof  equation  (2)). 
M  ore  specifically,  the  physics  of  the  photon  propagation  is  approximated  with  the  FEM  by 

A,,j  =  (-DVfj  ■  Vi pt  -  (4) 

where  fij  are  the  volume  elements  that  discretize  the  imaging  domain  and  form  a  geometrical 
mesh  defined  over  the  entire  imaging  domain.  A  is  integrated  over  this  imaging  domain.  The 
source  (in  this  case,  the  light  emitted  from  the  phosphors  which  were  excited  by  the  x-ray 
source)  and  boundary  integral  are  approximated  with  the  FEM  by 

bj  =  -  (  )  +<*y (f)  fjtjrids.  (5) 

\j—l  I  j—i  -/boundary 

The  time  component  of  the  luminescence  lifetime  is  ignored  since  the  measurements  in  this 
work  are  from  an  integrating  CCD  camera,  and  the  measurement  time  is  much  greater  than 
the  luminescence  lifetime;  effects  from  the  minimal  afterglow  of  the  phosphors  are  ignored. 

The  FEM  model,  G,  is  used  to  generate  estimates  for  the  photon  fluence  given  the  optical 
properties  of  the  tissue,  the  concentration  of  phosphors,  c,  and  the  FEM  mesh.  A  n  estimate 
for  the  photon  fluence,  </>,  can  be  calculated  by  solving  equation  (3): 

c/)  =  G  =  [A]-1{b}.  (6) 

In  this  paper,  the  imaging  domain  is  known,  and  is  assumed  that  the  endogenous  optical  tissue 
properties  are  known,  so  the  model  is  dependent  only  on  the  unknown,  c.  Figure  2  shows  the 
optical  photon  fluence  for  a  numerical  phantom  with  100:1  phosphor  concentration  between 
an  object  and  the  background.  Two  different  source  configurations  are  shown,  each  with  the 
x-ray  direction  of  propagation  in  the  horizontal  (left/right)  direction.  In  figures  2(a)  and  (b), 
the  x-ray  source,  indicated  by  the  red  circle,  irradiates  a  horizontal  line  passing  through  the 
background,  whereas  in  figures  2(c)  and  (d),  the  x-ray  irradiates  the  horizontal  line  passing 
through  the  middle  of  the  phosphor-containing  object. 
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FigureZ  (a)  Numerical  phantom  with  a  single  object  located  at  the  center  of  the  domain,  and 
a  10:1  contrast  in  phosphor  between  the  object  and  background.  The  x-ray  source  irradiates  the 
domain  along  the  long-axis  of  the  phantom  at  a  position  of  X=  10 ,Y=  -60:60,  (b)  The  emission 
fluence  field  of  (a),  (c)  Numerical  phantom  with  the  x-ray  source  irradiating  at  a  position  of  X  = 
40,  Y  =  -60:60,  (d)  The  emission  fluence  field  of  (c). 


2.2.3.  X-ray  nanophosphor  concentration  reconstruction.  The  goal  of  X  LT  is  to  determine 
the  phosphor  distribution.  This  spatial  distribution  can  be  determined  by  minimizing  the 
difference  between  the  measured  photon  flux  from  the  camera,  <t>M(r),  and  the  simulated 
photon  flux,  <ps[r)  at  identical  sample  locations.  This  is  accomplished  by  minimizing  the 
L2-norm  of  the  objective  function  in  an  optimization  routine: 

£2  =  (0s  -  <Pm)2  (7) 

where  £2  is  the  objective  function  to  minimize.  Because  of  the  large  dynamic  range,  the 
first  term  in  (7)  is  formulated  from  the  natural  logs  of  the  photon  fluxes.  This  is  an 
underdetermined  problem,  as  measurements  are  made  only  at  the  boundary.  Because  this 
problem  is  underdetermined,  the  model,  G,  is  linearized  with  a  Taylor  approximation  and 
formed  into  an  iterative  algorithm  as 

G(ci)  =  G(ci_i)+G,Ac  (8) 

where  G'  is  the  partial  differential  of  the  model  with  respect  to  the  concentration,  also  known 
as  the  J  acobian,  J .  M  inimizing  equation  (7)  with  respect  to  c  and  substituting  G(Cj)  from 
equation  (8)  into  <ps  and  J  for  G'  yields 

27  ((G(c,_i)  +  J  Ac)  —  4>m)  =  0 

Solving  for  the  concentration  yields 

Ac  =  —[JTJ]~lJT(<Ps  —  <Pm) 


(9) 
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This  problem  is  ill-posed,  so  it  is  solved  using  the  Levenberg-M  arquardt  ( 1963)  algorithm, 
which  includes  a  stabilization  parameter,  A,  in  the  inversion  to  avoid  singularities: 

A  c  =  -[JtJ+M]-1Jt(<Ps~4>m).  (10) 

Equation  (5)  is  iterated  until  a  minima  is  reached  (the  L2  norm  of  the  update  is  less  than  1% 
of  the  previous  iteration),  or  until  15  iterations  are  performed,  whichever  occurs  earlier.  The 
stabilization  parameter  is  reduced  at  each  iteration  as  the  algorithm  approaches  the  minimum 
and  converges  on  the  solution. 

Although  this  study  focused  on  applications  where  a  single  angle  is  ideal,  note  that 
this  algorithm  is  not  limited  to  a  single  angle.  Thus,  this  algorithm  is  appropriate  for  any 
sparse-angle  geometry. 

2.3.  Phantom  study 

The  performance  of  the  experimental  setup  and  the  reconstruction  algorithm  were  tested  by 
varying  the  size  and  location  of  a  lesion-simulating  object.  The  relationships  between  source- 
object  and  detector- object  distance  on  resolving  an  object  of  various  sizes  were  determined 
with  both  numerical  and  experimental  phantoms.  The  metrics  used  to  determine  system 
performance  were  object  location  and  object  size.  Location  error  in  both  the  x-ray  excitation 
and  optical  read-out  dimensions  was  determined  by  calculating  the  distance  between  the  true 
centroid  of  the  lesion  and  the  location  of  the  maximum  recovered  value  of  the  phantom. 
Object  size  in  both  the  x-ray  excitation  and  optical  read-out  dimensions  were  determined 
by  calculating  the  full  width  at  half  maximum  (FWHM  )  of  the  object  in  these  dimensions. 
Concentration  sensitivity  and  contrast  recovery  were  also  examined. 

2.3.1.  Numerical  phantoms.  Numerical  phantoms  wereutilized  to  testthe  position  accuracy, 
object-size  accuracy,  and  sensitivity  of  the  algorithms.  The  position  and  object-size  phantoms 
investigated  the  recovery  of  simulated  tumors  with  10:1  contrast  between  the  object  and 
the  background  and  phosphor  concentrations  of  10  /ug  ml-1.  These  lesions  varied  in  size 
between  2  and  14  mm,  and  were  placed  at  different  locations  in  the  phantom,  as  depicted  in 
figures  3(a)  and  4(a)  (note  the  multiple  arrows,  which  indicate  the  minimum  and  maximum 
extent  of  the  object  locations  investigated).  Sufficient  dose  (1  cGy)  was  given  to  yield 
signal-to-noise  (SNR)  greater  than  10  for  each  object  position— this  methodology  allowed  a 
performance  test  of  the  algorithm  for  all  object  positions. 

The  sensitivity  phantom  included  an  object  with  varying  concentration  (figure  5(a)),  and 
varying  contrast  (figure  5(b)).  The  6  mm  diameter  object  was  placed  at  the  center  of  the 
phantom  along  the  dimension  of  the  detectors  (the  long-axis),  and  moved  at  various  depths 
away  from  the  detectors.  We  used  the  phosphor  properties  from  K  andarakis  et  al  (1996)  to 
obtain  quantification  of  the  emitted  light  efficiency  for  their  lanthanum  oxysulfide:terbium 
phosphor,  which  was  1.39  x  1015  optical  photons/(Gy  x  mg).  We  incorporated  solid- 
angle  losses  as  well  as  losses  due  to  lens  inefficiency  (DO-1795,  Navitar  Imaging  Solutions, 
Rochester,  NY ).  SNR  below  10  was  assumed  to  be  too  low  to  detect. 

All  phantoms  were  two  dimensional,  and  measured  12  cm  x  6  cm.  Detectors  were  placed 
along  the  long-axis  of  the  phantom,  while  the  collimated  x-ray  source,  1  mm  wide,  was  scanned 
along  the  short  axis.  The  phosphor  used  for  this  experiment  mimicked  GOS:Eu,  demonstrated 
in  x-ray  luminescence  imaging  in  a  previous  study,  which  has  a  strong  luminescence  emission 
at  ~618  nm.  B  ackground  optical  properties  were  similar  to  that  of  breast  tissue  (Peters  et  al 
1990)  (na  =  0.0027,  At'  =  0.717). 
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Figure 3.  (a)  Numerical  phantom  experiment  examining  the  effect  of  moving  a  variable-sized 
object  away  (four  positions,  indicated  by  the  arrows)  from  the  x-ray  source  plane  (S),  while 
keeping  the  distance  from  the  detection  plane  (D)  fixed,  (b)  Recovered  object  location  error,  (c) 
Recovered  object  diameter  FWHM  with  respect  to  the  x-ray  source  dimension,  (d)  Recovered 
object  diameter  FWHM  with  respect  to  the  optical  detection  dimension. 

2.3.2.  Experimental  phantom.  The  experimental  phantom  is  shown  in  figure  6.  Figures  6(a) 
and  (b)  show  the  relative  layout  of  the  objects  in  the  phantom,  while  figure  6(c)  shows  the 
camera-eye  view  of  the  phantom.  Figure  6(d)  shows  an  image  of  the  phantom  while  the 
phantom  is  irradiated  by  the  x-ray  beam.  The  optically  clear  acrylic  phantom,  measuring 
12  cm  x  6  cm,  was  filled  with  India  ink  to  mimic  optical  absorption,  and  intralipid  to  mimic 
optical  scatter,  at  the  appropriate  concentrations.  The  optical  properties  were  p.a  =  0.0027  and 
Ai'  =  0.717,  as  determined  from  a  diffuse  optical  spectroscopy  system.  GOS:Eu  phosphor  at 
a  concentration  of  10  mg  mh1  was  added  to  two  cylindrical  inclusions,  one  5  mm  in  diameter, 
and  one  10  mm  in  diameter,  which  were  located  3  and  9  cm  from  the  edge  nearest  the  x-ray 
source,  respectively.  The  inclusions  were  both  imaged  at  various  depths  from  the  edge  nearest 
the  detector:  10,  15,  20,  30  mm;  these  dimensions  are  shown  more  clearly  in  figure  7(a). 
The  exposure  times  and  gains  were  1.5  s  at  gain  500,  3.5  s  at  gain  800,  3  s  at  gain  1000,  and 
7.5  s  at  gai  n  1000,  f or  i ncreasi  ng  depth.  D  ose  to  the  phantom  vari ed  dependi  ng  on  the  phosphor 
depth,  so  that  a  high  SN  R  could  be  acquired  while  the  phosphors  were  irradiated.  The  doses 
to  the  phantom  were  6.7, 15.6, 13.4,  and  33.4  cGy,  for  increasing  depth. 

3.  Results 

3.1.  Numerical  phantom  results 

As  described  above,  a  lesion-simulating  object  was  placed  in  various  locations  in  the  volume 
so  that  the  performance  of  LAX  LT  could  be  analyzed.  The  lesion  location  was  varied  with 


Limited-angle  x-ray  luminescence  tomography 


9 


Detector  to  Object  Distance  (mm) 


— e- 

-  2mm 

- B— 

-  4mm 

— 0- 

-  6mm 

— v- 

-  10mm 

14mm 

Figure4  (a)  Numerical  phantom  experiment  examining  the  effect  of  moving  a  variable-sized 
object  away  (eight  positions,  indicated  by  the  arrows)  from  the  optical  detection  plane  (D),  while 
keeping  the  distance  from  the  x-ray  source  plane  (S)  fixed,  (b)  Recovered  object  location  error, 
(c)  Recovered  object  diameter  FWHM  with  respect  to  thex-ray  source  dimension,  (d)  Recovered 
object  diameter  FWHM  with  respect  to  the  optical  detection  dimension. 
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-  b  -  Truth:  10:1 
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Figured  (a)  Numerical  phantom  experiment  examining  the  effect  of:  (a)  variable  object  to 
background  contrast  and  (b)  variable  concentration,  of  a  6  mm  diameter  object  versus  depth  from 
the  detection  plane. 


respect  to  the  source-axis  and  detection-axis  separately  to  determine  the  effects  of  source- 
object  distance  and  detector- object  distance  on  the  ability  to  resolve  the  object.  Figure  3 
shows  the  results  of  maintaining  a  fixed  detector- object  distance  while  varying  the  source- 
object  distance  (varying  depth  with  respect  to  the  x-ray).  These  dimensions  are  depicted  in 
figure  3(a).  The  location  error  for  the  2-14  mm  objects  is  shown  in  figure  3(b).  This  result 
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Figure 6  Experimental  phantom:  (a)  and  (b)  two  views  from  above  the  black-box  including 
the  phantom  and  collimator  bricks,  (c)  Image  taken  by  the  high-sensitivity  CCD  camera  of  the 
phantom  under  ambient  light  and  no  irradiation,  (d)  CCD  image  with  no  ambient  light  and  x-ray 
irradiation. 


demonstrates  the  excellent  ability  of  the  algorithm  to  spatially  resolve  an  object,  as  all  errors 
in  distance  are  lower  than  3.5  mm.  Figure  3(c)  shows  the  recovered  object  size  in  the  direction 
of  the  x-ray  excitation,  which  demonstrates  the  insensitivity  with  respect  to  depth,  with  an 
average  standard  deviation  in  the  error  in  recovered  size  of  0.82  mm  for  varying  source-object 
distances.  In  figure  3(d),  the  recovered  size  of  the  object  in  the  optical  dimension  is  examined 
as  a  function  of  source-object  depth.  The  algorithm  is  able  to  distinguish  the  varying  sizes  of 
the  objects,  with  a  standard  deviation  in  the  error  in  recovered  size  of  0.83  mm  for  the  varying 
source-object  distance.  There  is  a  slight  tendency  for  blurring  in  the  optical  dimension  due  to 
the  scattering  of  the  optical  photons. 

Figure  4  shows  the  effects  of  varying  the  depth  of  an  object  with  respect  to  the  optical 
dimension.  Figure  4(a)  shows  the  physical  dimensions  of  the  phantom.  The  object  was  moved 
from  0.5  to  4.5  cm  from  the  detector,  while  remaining  at  a  fixed  distance  from  the  x-ray  source. 
Figure  4(b)  shows  the  recovered  location  error  of  the  centroid  of  the  object  with  respect  to 
depth  from  the  optical  detector.  In  this  experiment,  the  ability  to  resolve  an  object  of  ~2  mm 
was  at  the  limit  of  the  system  with  the  1  mm  collimated  x-ray  beam  that  was  used.  Flowever, 
the  objects  sized  4  mm  and  larger  were  resolved  with  higher  accuracy.  As  expected,  location 
error  increases  slightly  as  the  depth  increases,  and  larger  objects  are  resolved  more  accurately 
than  smaller  objects.  The  advantage  of  using  the  x-ray  is  highlighted  here,  as  even  at  4.5  cm  in 
depth,  the  object  is  resolved  with  less  than  1  cm  total  error.  Figure  4(c)  shows  the  recovery  of 
object  size  in  the  x-ray  dimension,  similar  to  figure  3(c).  The  mean  standard  deviation  for  the 
error  in  recovered  object  size  was  1.36  mm.  Figure  4(d)  shows  that  the  ability  to  determine 
the  size  of  the  object  in  the  optical  dimension  has  increased  variability  compared  to  the  x-ray 
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Figure 7.  (a)  Experimental  phantom  of  two  objects  moved  away  from  the  optical  detection  plane 
(D)  at  four  positions,  (b)  Recovered  object  location  error,  (c)  Recovered  object  diameter  FWHM 
with  respect  to  the  x-ray  source  dimension,  (d)  Recovered  object  diameter  FWHM  with  respect  to 
the  optical  detection  dimension. 


dimension,  as  the  mean  standard  deviation  of  the  error  in  recovered  size  with  respect  to  depth 
is  2.4  mm.  This  larger  standard  deviation  is  consistent  with  object  blurring  at  depth  from  the 
optical  detector;  this  is  expected  due  to  the  ill-posed  nature  of  the  algorithms  needed  for  the 
optical  photon  modeling.  Similar  to  figure  3(d),  this  system  is  able  to  distinguish  between 
the  different- si  zed  objects  in  the  optical  dimension,  although  again,  there  is  a  tendency  for 
dilation  in  the  optically  resolved  dimension. 

Overall,  the  numerical  phantom  results  show  localization  accuracy  with  median  error  of 
2.2  mm  (mean  of  6.3  mm),  or  4.1%  (mean  of  11.5%)  of  object  depth  for  all  lesions.  The 
recovered  lesion  size  has  lower  size-bias,  with  median  error  of  -8.1%  versus  87.5%  (mean  of 
19.4%  and  118.3%,  respectively),  inthex-ray  excitation  direction  versus  the  optical  direction, 
respectively.  Again,  this  optical  dilation  is  expected  due  to  the  increased  optical  scatter 
compared  to  x-ray.  This  technique  is  invariant  in  recovered  size  with  respect  to  depth,  as  the 
standard  deviation  is  less  than  2.5  mm. 

The  concentration  phantom  is  shown  in  figure  5(a).  The  required  dose  (in  Gy)  to  reach  an 
SNR  of  10  is  plotted  for  varying  concentrations.  It  is  apparent  from  this  calculation  that  a  6 
mm  diameter  object  with  ytg  mb1  concentration  is  detectable  in  this  geometry  with  standard 
CT  doses.  With  the  doses  currently  used  in  IORT,  concentrations  to  ng  ml-1  are  detectable  at 
depth.  As  shown  in  figure  5(b),  with  sufficient  signal,  contrast  can  be  recovered  for  all  contrast 
to  background  ratios  tested,  for  depths  up  to  45  mm.  Here,  the  advantage  of  the  ability  of  the 
collimated  x-ray  to  selectively  excite  the  phosphors  is  clear. 


3.2.  Experimental  phantom  results 

Photographs  of  the  experimental  phantom  are  shown  in  figure  6;  the  schematic  of  the 
experimental  phantom  is  shown  in  figure  7(a),  and  the  results  are  shown  in  figures  7(b)- (d). 
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Figures  Image  reconstructions  of  two  objects  with  depths  of  (a)  10,  (b)  15,  (c)  20  and  (d)  30  mm 
from  the  detection  plane. 


Reconstructions  for  each  case  are  shown  in  figure  8,  with  increasing  depth  with  each  row.  In 
this  experiment,  a  1-5  mm  object,  and  a  1-10  mm  object  were  imaged  at  various  depths.  For 
all  depths,  the  recovered  location  error  for  both  objects  is  less  than  6  mm,  and  is  independent 
of  object  size.  Positional  errors  are  recovered  within  8.6%  of  the  effective  depth  for  a  5  mm 
object,  and  within  5.2%  of  the  depth  for  a  10  mm  object.  Similar  to  the  results  from  the 
numerical  phantom,  the  ability  to  discern  object  size  is  highly  accurate  in  the  x-ray  dimension, 
and  is  invariant  with  depth.  The  ability  to  resolve  the  object  in  the  optical  dimension  is  accurate 
to  within  2  mm.  Object-size  median  error  is  within  2.3%  and  2%  for  the  5  mm  and  10  mm 
objects,  respectively. 
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Figure 9.  (a)  Object  location  linearity  of  the  depth  of  the  object  with  respect  to  the  detection 
plane,  (b)  Object  location  linearity  of  both  the  depth  of  the  object  with  respect  to  the  detection 
plane  summed  with  distance  with  respect  to  the  source  plane. 


4  Discussion 

Molecular  imaging  has  long  been  identified  as  potentially  a  vital  tool  in  surgery  to  aid  in 
the  identification  of  important  tissue  structures  and  tumor  tissue  (Gregorie  etal  1968).  U  ntil 
recently,  the  use  of  molecular  imaging  in  surgery  has  been  limited,  both  due  to  the  lack  of 
appropriate  imaging  tools,  as  well  as  the  lack  of  highly  specific  contrast  agents  (the  problem 
with  most  endogenous  fluorescence  agents  in  the  body).  Recent  developments  in  both  of  these 
areas  have  increased  the  interest  of  molecular  surgical  guidance  (Weissleder  and  Pittet  2008). 
For  instance,  fluorescence  molecular  tomography  (FMT)  shows  great  promise  in  positively 
affecting  surgical  outcomes,  as  it  is  non-ionizing,  has  a  wealth  of  knowledge  of  contrast 
agents,  and  can  be  implemented  at  low  cost.  FMT  has  found  utility  in  tumor-margin  excision 
(Stummer  etal  2000),  sentinel  lymph-node  mapping  (Tanaka  etal  2006),  and  avoidance  of 
critical  structures  such  as  nerves  (Boyette  etal  2007).  Alternatively,  sentinel  lymph-node 
identification  with  radio-emission  imaging  is  widely  adopted  (Alex  and  Krag  1993);  its  use  in 
surgical  resection  is  not  well  established. 

The  lack  of  wide  adoption  of  molecular  imaging  in  surgical  applications  can  partly  be 
blamed  on  the  limited  flexibility  of  these  imaging  technologies.  While  FMT  is  ideal  for 
structures  at  the  surface  and  several  millimeters  beyond,  its  limitations  in  depth  penetration 
preclude  its  ability  to  unmask  disease  that  exists  several  centimeters  beneath  the  surface. 
Radio-emission  imaging  is  limited  because  of  the  need  to  image  gamma  photons  which  are 
highly  penetrating  in  tissue— this  prevents  the  identification  of  lesion  depth,  and  increases  the 
burden  on  the  surgeon  to  patiently  dig  through  the  tissue  until  the  object  is  identified.  In  the 
case  of  a  sub-millimeter-sized  tumor  tissue,  this  practice  would  be  overly  laborious. 

Clearly,  a  void  exits  in  the  ability  of  molecular  instrumentation  to  resolve  millimeter  or 
sub-millimeter  objects  at  depthsof  several  centimeters.  LAXLT  may  be  able  to  fi  1 1  this  voi  d.  Its 
advantage  is  the  ability  to  resolve  objects  at  several  centi  meters  of  deep.  F  igure  9  demonstrates 
this  ability,  as  it  is  shown  that  both  the  5  mm  and  10  mm  objects  in  the  experimental  phantom 
were  able  to  be  resolved  at  a  detector- object  depth  of  beyond  3  cm.  This  is  a  stark  contrast 
to  FMT,  as  Kepshire  etal  (2007)  demonstrated  that  depth  linearity  degrades  beyond  1  cm.  If 
the  depth  penetration  of  the  x-ray  source  is  considered,  LAX  LT  is  linear  beyond  10  cm.  This 
increased  performance  is  due  to  the  low  scatter  and  high  penetrability  of  the  x-ray  excitation, 
which  allows  more  ideal  imaging  geometries  to  be  chosen  for  surgical  guidance. 
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The  benefit  of  the  technique  developed  in  this  study  compared  to  previous  developments 
with  X  LT  is  threefold:  increased  temporal  resolution  is  possible  because  of  the  optical  read¬ 
out  of  two  dimensions,  which  eliminates  x-ray  encoding  in  those  dimensions:  lower  dose  is 
possible  because  of  the  decreased  x-ray  excitation.  M  ost  importantly,  this  technique  is  suitable 
for  geometries  where  full  angular  x-ray  encoding  is  not  possible.  Therefore,  this  technique 
enables  X  LT  to  be  used  for  surgeries  such  as  breast  lumpectomy:  we  foresee  applications  such 
as  this  to  be  one  of  the  most  important  future  applications  in  imaging.  In  breast  lumpectomy,  a 
full-angular  encoding  with  x-ray  is  not  desirable  because  deep  critical  structures  are  irradiated. 
By  instead  implementing  a  limited-angle  technique,  only  the  breast  may  be  irradiated.  For 
lumpectomy,  this  technique  would  be  desirable  to  verify  the  position  of  the  lesions  in  the 
surgical  supine  position,  and  to  visualize  remaining  disease  after  resection. 

In  comparing  this  technique  to  other  image-guided  surgical  modalities  such  as  FM  T,  the 
advantages  of  LAX  LT  are  clear.  LAXLT  was  able  to  resolve  4  mm  objects  at  detector- object 
depths  greater  than  4.5  cm,  compared  to  FM  T  which  is  limited  in  to  about  1  cm  in  resolving 
object  dimensions  (K  epshi  re  eta/  2007).  The  most  impressive  aspect  indicated  by  these  results 
is  the  reduction  in  the  blurring  of  the  object  with  respect  to  depth:  this  highlights  the  advantage 
of  incorporating  the  x-ray  excitation,  which  due  to  its  relatively  insignificant  scatter  at  these 
depths  can  pinpoint  the  depth  of  the  object  with  respect  to  the  optical  depth  dimension.  The 
optical  read-out  is  then  used  to  determine  the  other  dimensions.  Since  the  depth  is  known  to 
high  accuracy  (shown  in  figure  4),  the  diffuse  algorithm  properly  models  the  diffuse  nature 
of  the  light,  and  significantly  reduces  object  blurring,  as  presented  in  figure  5.  Thus,  with 
LAX  LT,  dose  and  imaging  time  are  reduced  significantly  compared  to  X  LT. 

It  is  intuitive  that  the  resolution  will  be  limited  by  the  width  of  the  x-ray  beam,  at  shallower 
depths.  At  deeper  depths,  the  scatter  of  the  x-ray  beam  should  be  taken  into  account.  In  this 
study,  a  1  mm  collimated  x-ray  beam  was  able  to  successfully  visualize  a  2  mm  object  up  to 
2  cm,  and  a  4  mm  object  up  to  4.5  cm.  H  igher  resolution  should  be  attained  with  a  narrower 
beam:  we  are  currently  investigating  this  effect.  The  ability  to  resolve  an  object  at  depth  is 
instead  limited  by  dose. 

A  Ithough  this  study  was  adapted  to  geometries  best  suitable  to  tumor  resection,  it  should 
be  noted  that  the  algorithm  presented  in  this  work  is  generalizable  to  any  geometry:  it  is 
especially  useful  for  sparse-angle  geometries.  An  additional  use  for  this  technique  would 
be  for  intraoperative  probes,  where  simple  coregistration  (assuming  the  catheter  is  radio¬ 
opaque)  between  the  x-ray  source  and  the  optical  detection  catheter  could  provide  assessment 
of  molecular  status  at  a  remote  location. 

An  increase  in  dose  to  potentially  healthy  tissue  is  one  disadvantage  of  this  technique. 
Although  we  calculated  that  the  maximum  dose  to  the  phantom  was  33.4  cGy,  this  dose  will 
be  reduced  to  6  cGy  with  a  more  favorable  optical  setup  where  the  camera  is  closer  to  the 
object.  This  dose  may  be  reduced  further  with  a  contact  optical  setup.  Still,  this  dose  may 
preclude  its  use  in  a  screening  setting,  especially  due  to  the  increasing  awareness  of  increasing 
radiation  exposure  in  medicine  (Caoili  etal  2009).  Flow  ever,  post-surgical  radiation  therapy 
is  commonly  prescribed  as  a  means  to  destroy  cancer  cells  that  may  not  have  been  removed 
around  the  margin  (Dirbas  2009,  M  unshi  2007).  In  accelerated  partial  breast  irradiation, 
5-20  Gy  of  radiation  is  given  to  the  resected  cavity  in  1-5  fractions,  to  reduce  the  morbidity 
of  whole  breast  radiation  therapy  (Ross  2005).  In  this  context,  LAXLT  may  have  great 
utility  to  identify  larger  regions  that  may  have  been  missed  during  surgery,  and  may  require 
subsequent  surgical  investigation.  Its  sensitivity  at  this  dose  might  enable  the  identification  of 
micro-disease,  or  a  very  low  concentration  to  be  injected. 

The  other  disadvantage  of  this  technique  is  its  use  of  nanoparticles.  Although  we  have 
demonstrated  low  toxicity  of  our  nanoparticles  in  cells,  nanoparticles  will  have  different 
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effects  in  a  living  system.  This  topic  is  beyond  the  scope  of  this  paper,  but  it  is  important  to 
recognize  that  cancer  nanotechnology  is  a  major  venture  in  the  National  Cancer  Institute,  and 
this  technique  will  benefit  from  the  knowledge  gained  from  this  program. 


5.  Conclusions 

LAXLT  has  been  developed  as  a  means  to  image  molecular  deeper  than  is  available  with 
FMT.  A  reconstruction  algorithm,  based  on  a  hybrid  x-ray  excitation/diffuse  optical  emission 
model,  was  tested  in  a  numerical  and  experimental  phantom  that  had  dimensions  similar  to 
the  human  breast.  It  was  found  that  objects  as  small  as  2  mm  in  diameter  could  be  resolved 
at  depths  of  up  to  4-5  cm.  It  was  shown  that  the  depth  of  the  object  with  respect  to  the  x-ray 
source  position  had  little  effect  on  object  recovery  in  this  volume.  It  was  then  demonstrated 
experimentally  that  both  a  5  and  10  mm  object  could  be  resolved  at  depths  of  at  least  3  cm; 
these  results  agreed  with  the  numerical  phantom  results,  thus  validating  the  simulation.  If  the 
challenges  to  engineering  biocompatible  phosphors  can  be  resolved,  LAXLT  may  have  utility 
in  surgical  applications  where  small  lesions  must  be  imaged  at  depth  of  a  few  centimeters, 
such  as  during  breast  lumpectomy  surgery. 
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